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We introduce an exactly solvable SU(2)-invariant spin-1/2 model with exotic spin excitations. 
With time reversal symmetry (TRS), the ground state is a spin liquid with gapless or gapped spin-1 
but fermionic excitations. When TRS is broken, the resulting spin liquid exhibits deconfined vortex 
excitations which carry spin-1/2 and obey non-Abelian statistics. We show that this SU(2) invariant 
non-Abelian spin liquid exhibits spin quantum Hall effect with quantized spin Hall conductivity 
o~xy — h/2ir, and that the spin response is effectively described by the SO(3) level-1 Chern-Simons 
theory at low energy. We further propose that a SU(2) level-2 Chern-Simons theory is the effective 
field theory describing the topological structure of the non-Abelian SU(2) invariant spin liquid. 



Quantum spin liquids are zero temperature exotic 
states of magnets which exhibit no classical order. Their 
existence is often due to strong quantum fluctuations 
and/or geometric frustrations. In two dimensions (2D), 
the notion of quantum spin liquids were first proposed by 
Anderson in 1973 as the candidate ground state of the 
spin-1/2 AFM Heisenberg model on the triangular lattice 
In addition, soon after the discovery of the cuprate 
high temperature superconductors quantum spin liquids 
are proposed as the root states of the high T c supercon- 
ducting states [2 4]. However, up to now there seems no 
definitive proof of the existence of real spin liquid mate- 
rials for D > 2 even though a few promising candidates 
have emerged recently Consequently, further ex- 

perimental scrutinies are needed to firmly establish the 
existence of real spin liquids in D > 2. Of course, mi- 
croscopic theoretical models with reliable solution of spin 
liquid ground states are also desired. 

In order for the elementary excitations to carry well- 
defined spin quantum number, SU(2) invariance is im- 
portant. Recently, a number of exactly solvable SU(2) 
invariant models [7H12J with spin liquid ground states 
have been proposed. They have the common trait that 
the low energy /long wave length physics is captured by 
the short-range resonating valence bond (RVB) @, 

In this paper, we take a different route to construct 
a SU(2) invariant model with spin liquid ground states 
by generalizing the recently discovered Kitaev model (l3l | 
in a SU(2) symmetric way. Our model [see Eq. (JTJ] is 
a spin-1/2 system on the decorated honeycomb lattice 
(also known as the star or 3-12 lattice, see Fig. [1]). Our 
spin model can be reduced to three species of free Majo- 
rana fermions coupled to background Z 2 gauge field such 
that the model is exactly solvable, in similar spirit of 
the original Kitaev model that have only one species of 
free Majorana fermions. In the presence of time reversal 
symmetry (TRS), it exhibits a spin liquid ground state 
with gapless or gapped spin-1 excitations, depending on 
parameters in the model. Strikingly, these spin-1 exci- 
tations obey fermionic statistics. When TRS is broken 



cither spontaneously or explicitly, each flavor of Majo- 
rana fermions behave like the Bogoliubov quasiparticle 
of a chiral p + ip superconductor. As a result, a Z 2 vor- 
tex binds three Majorana zero modes, protected by the 
SU(2) symmetry. Due to the existence of the odd num- 
ber of Majarana zero modes, the vortex excitations obey 
non-Abelian statistics, which might find potential appli- 
cation in topological quantum computing [141 ] . We fur- 
ther show that the non-Abelian vortex excitations carry 
spin-1/2 quantum number. As far as we know, this is 
the first realization of non-Abelian spinon in an exactly 
solvable SU(2) symmetric spin model. Note that a SU(2) 
symmetric wave function with non-Abelian spinon excita- 
tions was proposed recently in a spin-1 system[15[ but it 
is not clear which Hamiltonian can realize it as a ground 
state. 

In the SU(2) invariant non-Abelian phase, the system 
exhibits a spin quantum Hall effect [161 ] with quantized 
spin Hall conductivity a^, y = H/2n (which is twice of that 
in a d + id superconductor (it|.) The quantized spin Hall 
response is shown to be described by the SO (3) level-1 
Chern-Simons gauge theory at low energy. Furthermore, 
according to the threefold ground state degeneracy, non- 
Abelian statistics of vortices, and the SU(2)-invariancc 
of the ground states, we propose that the low energy 
topological field theory for the non-Abelian phase is the 
SU(2) level-2 Chern-Simons theory 0G3]- 




FIG. 1. The schematic representation of the decorated hon- 
eycomb lattice. The triangles are labeled by i,j, the sites 
within each triangle are labeled by 1,2,3, and the type of 
inter-triangle bonds is x, y, or z. 



2 



The model: We consider the following SU(2) invari- 
ant Hamiltonian on the lattice shown in Fig. [if a): 



H 



(1) 



A-link (ij) 



where i, j label the triangles and Si = S^i + S^2 + Si^ 
is the total spin of the zth triangle (Si >a is the spin-1/2 
operator on site a = 1, 2, 3 of the ith triangle.) The op- 
erators are defined as follows: rf = 2(Si i ■ Si 2 + 1/4), 
rf = 2(S M ■ S li3 - S,, 2 ■ S <|3 )/A and r? = 4S M ■ (S;, 2 x 
Si,3)/V3- The parameter J is the strength of the intra- 
triangle spin exchange coupling while J\ (A = x,y,z) de- 
scribes the inter-triangle couplings on the type-A links. 
Because [Sf.Sj] = and [Sf,r/] = 0, the operator Sf 
commutes with the Hamiltonian for all i [2(|. As a result, 
the total spin of each triangle is a good quantum number, 
so we can use them to subdivide the Hilbert space. 

For each triangle, three spin-l/2's can be decomposed 
in terms of their total spins: = §©l©5- Within 

the space spanned by \ © \ it is straightforward to check 
that n satisfy the SU(2) algebra [rf ,rf] = 2ie a ^rf as 
well as the Clifford algebra {rf ,rf} = 2S a ^ , which im- 
plies that r* are the Pauli matrices. In fact the four 
states in i © i can be labeled by \af = ±1, rf = ±1) with 
<?, = 2Si since [rf,Sf] = 0. Consequently we may view 
as pesudo-spin 1/2 which distinguishes the "orbital" 
degree of freedom within the two degenerate spin-1/2 
multiplcts. In the remaining spin-3/2 space, rf = rf = 
because the three spins are totally polarized for Si = 3/2; 
r 4 A do not satisfy SU(2) algebra. 

When J ^> J\, the ground state and the low lying 
excited states all lie in the sub-Hilbert space £0 where 
Si = 1/2 for all triangles. [It can be shown that the 
ground state must lie in £ as long as J > |( J x + J y + J z ).\ 
For simplicity, we will assume J J\ hereafter which 
allows us to focus on Cq. In this sub-Hilbert space, apart 
from a constant, the Hamiltonian in Eq. ([I} reduces to 



H = 



1 



A-l 



E 

ink (ij) 



(2) 



which seems complicated but is actually exactly solvable 
as shown below. To solve the model, we introduce Majo- 
rana fermions representations [l3j for the Pauli matrices 
erf and rf as follows: 



e —i4cl rt = - e -—id^dl3) 



where a, j3 = x,y,z and cf, df are Major ana fermions. 
As usual, the Majorana fcrmion representation is over- 
complete and the following constraint is needed for a 
physical wave function: |^ , ) p h ys i s m physical Hilbert 
space iff 



A|*> phys = |*) phV8 ,Vi 



phys : 



(4) 



where Di = —ic^c^c^d^d^d^. In other words, any state 



n, l 1 - 



+ D,- 



acted by the projection operator P 
physical state. 

In terms of Majorana fermion operators, it is straight- 
forward to rewrite the spin Hamiltonian Eq. ^) as follows 



H 



E JijUij [icfc* + ic\c y - + ir,, : j 

(ij) 



(5) 



, 3 - -id£d$ 
x,y,z) link (ij). 



where v 
(A 

cause [uijjH] = 0, and 



and Ji 



Ja/4 on the type-A 



It is clear that H = PUP. Be- 
[uij,Uiij>] = 0, Uij are good quan- 



tum numbers with eigenvalues ±1. It is obvious that 
Eq. ([5]) is invariant under the following local Z 2 gauge 
transformation cf — > Ajcf and Uij — > AiUijAj, Aj = ±1. 
Eq. ([5]) describes three species of free Majorana fermions 
coupled with background Z 2 gauge field Uij. In addition 
to the Z 2 gauge symmetry Eq. ([5]) has a global SO (3) 
symmetry which rotate among the three species of Ma- 
jorana fermions, which is the consequence of the SU(2) 
symmetry of the original spin model. 

Lieb's theorem [2l| requires that the ground state lies 
in the zero flux sector, namely <f> p = for every hexagon 



plaquette p where exp(i</> p ) = F 



{jk)e P lu 3k 



The zero flux 



sector is realized by choosing Uy = 1 with € A(B) 
sublattice. In the zero flux sector, it is straightforward 
to show that there are fermionic excitations [22j that are 
gapped or gapless with Dirac-like dispersion, as discussed 
below. Besides these fermionic excitations, vortex exci- 
tation on plaquette p is created when <p p = n. 

Gapless/gapped spin liquid: The spectrum 
of each species of Majorana fermions in Eq. §5§ 
can be obtained by a Fourier transform cf = 



E 



4>a. 



where 



A/B is the sublattice index of site i, k a ,il>p,u.',b} = 
SapSk kiSab, and the summation is over one-half of the 
Brillouin zone. The three species have identical dis- 
persion Ek = ±| E Q J a e lU ' ea I , which is gapless with 
Dirac-like dispersion for J x + J y > J z and is gapped 
otherwise. (Here e a label nearest-neighbor vectors.) It 
is straightforward to calculate the spin-spin correlation 
with (SfSj) ~ |rj — ij| -4 for |r, — r,-| » 1 in the gap- 
less phase but in exponential decay for the gapped phase. 
This is in contrast with the original Kitaev model and its 
variants studied previously where the spins are uncorre- 
cted beyond nearest neighbors 23, 2~3 |. 

Fermionic spin-1 excitations: In the zero flux sec- 
tor, what is the nature of fermionic excitations? We show 
below that these fermionic excitations carry spin-1 quan- 
tum number. For instance, by introducing the complex 
fermion operators /i j2 = (cf — icf)/2, the Hamiltonian 
Eq. ([5]) can be rewritten as 



H 



(ij) 



13 "ij 



(ifUi,.-ifl,fi,*) ( 6 ) 
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which indicates that fi tZ are free complex fermions and 
their fermion number is conserved. The dispersion of 
f z complex fermions and c z Majorana fermions is again 
£ k = ±| J^a ^e lke °| but k is over the full (half) Bril- 
louin zone for f z (c z ). Since 5? - crf/2 = fj j^ - 1/2, 

it is clear that the fermions created by fj z carry S z = 1, 
which can be loosely called fermionic magnons with 
S z = 1. The c z quasi-particles carry S z = 0. 

In the gapped phase, due to total f z fermion num- 
ber conservation (a result of total S z conservation), the 
Chern number of f z can be defined and it is zero. Simi- 
larly, a spectral Chern number can also be defined for c z 
Majorana fermions and it is also zero. Due to the zero 
Chern number, a vortex excitation has no fermion zero 
mode for f z and no Majorana zero mode for c z ; its spin 
quantum number is zero. 

Non-Abelian spinons: A spinon carries spin-1/2 
quantum number and is essentially half of a spin-1 ex- 
citation. Deconfincd spinon is commonly regarded as 
a hallmark of a quantum spin liquid. A spinon excita- 
tion with S z — 1/2 would carry half fermion number of 
f z . Excitations with half fermion quantum number are 
known to exist as fermion zero modes in some sort of de- 
fect such as domain walls and vortices. Thus, we study 
vortex excitations. Since a fermion zero mode with half 
fermion number is expected in a vortex with flux tt when 
the insulator carries Chern number ±1. 

To obtain a finite Chern number for f z , following Ref. 
[25j . we consider additional interactions which breaks 
TRS externally (TRS can b e sp ontaneously broken by 
"decorating" the star lattice |26(): 



H '=h y: 



T i T j T k 



[Si ■ S fc ] , (7) 



where (ij) a (jk}^ labels three neighboring sites, which 
are ordered in a clockwise way within the corre- 
sponding hexagon plaqucttc, forming a triad whose 
links are a and /3 respectively. In term of fermions, 



H ' = l E(y)o-fc) UijUjk (icfc x k + icfc y k + ic z c z k ) 




2% i 



I Efe>0fc) UijUjk[2(iftjk, z - ifl,Ji,z) + This 
next-nearest-neighbor hopping term generates a mass 
term that gaps out the Dirac point of f z fermions and 
the resulting insulator of f z has Chern number v = 1, 
which can be computed directly and is also indicated by 
the one gapless chiral edge mode of f z on a system with 
boundary, as shown in Fig. [2] Since f z fermions carry 
S z = 1, the finite Chern number v = 1 implies the spin 
quantum Hall effect J z ' 1 = a?jdjB z , where J z ' t is the cur- 
rent of spin polarized along the z-direction flowing along 
the i = x, y-dircction, with quantized spin Hall conduc- 
tivity <i s xy = = H/2tt. 

Because of the finite spin Hall conductivity, by insert- 
ing a 7r (or equivalently — tt) flux locally on plaquette 
p, namely creating a vortex excitation, the spin is cumu- 
lated around the vortex by the amount S z = a xy ■ (±tt) = 
±ft/2, implying that the vortex carries spin-1/2 quantum 
number. The sign ambiguity reflects the facts that there 
is one zero mode associated with the vortex excitation 
and that occupying it or not gives rise to S z = ±1/2. 
The spin-1/2 nature of a vortex excitation is further ver- 
ified numerically. For instance, by creating two well sep- 
arated vortices in a finite lattice and occupying the zero 
mode associated with each vortex, we obtain S z = 1/2 
for each vortex, as shown in Fig. [U(a). 

Similarly, the c z Majorona fermions have spectral 
Chern number ±1 and the edge has a gapless chiral Ma- 
jorana mode. Moreover, there is one c z Majorana zero 
mode associated with each vortex excitation [13j . Due to 
this unpaired Majorana zero mode from c z , the vortex 
excitations obey non-Abelian statistics. In the following 
we refer to them as "non-Abelian spinons" since they also 
carry spin-1/2 quantum number. It is now clear that a lo- 
cal vortex excitation actually binds three Majorana zero 
modes, which arc treated "on bias" as a complex fermion 
zero mode plus a Majorana zero mode above. The three 
Majorana zero modes will not mix and split because of 
the SU(2) symmetry. 

We further computed the total energy of a pair of vor- 
tices as a function of the inter- vortex distance and the 
result is shown in Fig. Elb). The fact that the energy 
decreases with the increasing distance indicates that the 



S<(v,,) 




two spinon energy 




(a) 



(b) 



FIG. 2. The energy spectrum as a function of k y for f z 
fermions on a cylinder periodic in y-direction. There is ex- 
actly one gapless chiral edge state. The parameters used for 
both Fig. 2 and Fig. 3 are J x — J y — J z = 1 and h = 0.3. 



FIG. 3. (a) Two localized spinons each with S z — 1/2 are 
realized by creating two vortex excitations, (b) The energy 
of creating two localized spinons as a function of distance x 
between then. 
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vortex-vortex interaction is repulsive and creating two 
far separated spinons only cost a finite energy. The non- 
Abelian spinons are deconfined. 

Topological field theory: A topological phase is 
generically described by a topological field theory In 
the non-Abelian spin liquid discussed above, we expect 
that the finite spin quantum Hall effect have a topological 
field theory description which we will derive below. For 
simplicity, we assume J x = J y = J z hereafter. Without 
the h term, the dispersion is gapless with the Dirac point 
at K. (Note that each species of Majorana fermions has 
only one Dirac cone due to their Majorana nature.) A 
finite h term acts like a mass term for the Dirac fermions 
at K. In the continuum limit, the low energy physics is 
then described by the following Euclidean action 



S = dtd 2 x^\i^d i , + imjilj, 



(8) 



where ip = (ip x A,^xB,ip y A,'>pyB,'>p z A,tpzB) T and ip = 
iip'^f . Here m is the mass and j > x ' y = I®a z ' v ' x where I 
is a 3 by 3 identity matrix with vector indices x, y, z and 
a Pauli matrices with sublattice indices. 

It is clear that the three species of massive Dirac 
fermions possess a global SO(3) symmetry, which is in- 
herited from the SU(2) symmetry of the original spin 
model. The continuous SO (3) symmetry allows us to in- 
troduce external spin gauge fields , which couple with 
the spin current J a '' = ipH' 1 ^ 1 *'!!;, where [t a ]b c = ie abc , in 
the following way: 



dtd 2 xip(t,x.) i^(d^ + iA a t a ) + im ^(i,x).(9) 



By integrating out fermions, we obtain an induced action 
for A^ , whose lowest-order imaginary part is a topological 
term: 



S eS [A] = ^- I dtd 2 xe^ x 



1 



automatically. Without derivations, we pro- 
pose that the following SU(2) lcvcl-2 [SU(2) 2 ] Chern- 
Simons theory 0, 



i— l dtd 2 x^ vX 
4ir 



(12) 



is the low-energy effective theory for the non-Abelian 
spin liquid phase. There are several reasons for such a 
proposal. First it is natural to use SU(2) as the gauge 
group; the Majorana fermions has both SO (3) and 
symmetries, which combine leading to an effective SU(2) 
as expected from the spin SU(2) symmetry of the model. 
Secondly, as discussed in Ref. [19| the edge theory corre- 
sponds to the SU(2)2 Chcrn-Simons theory [Eq. (|12p] is 
the chiral sector of SU(2)2 Wess-Zumino-Witten model 
with central charge c = 3/2, which is consistent with the 
fact that the edge theory here has three copies of chiral 
Majorana modes. Thirdly, the SU(2)2 Chern-Simons the- 
ory on a torus is threefold degenerate, which is identical 
with the degeneracy computed from our lattice model. 

Concluding remarks: We have shown that the ex- 
actly solvable SU(2)-invariant spin-I/2 model on the dec- 
orated honeycomb lattice exhibits quantum spin liquid 
ground states with fermionic magnons or non-Abelian 
spinons. Interestingly, a recently discovered material, 
called Iron Acetate [3l|, realizes a spin model on the 
decorated honeycomb lattice, which adds some hope that 
our model may be realized in similar family of materials. 
Moreover, we believe that the model can be potentially 
realized by loading cold atoms in specially designed op- 
tical lattices under appropriate circumstances [30(. 
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DE-AC02-05CHI123I. 



(10) 



which is a SO(3) level-I Chern-Simons action. This topo- 
logical term describes the spin responses of the system to 
the external spin gauge fields. Physically, A§ = B a is the 
external magnetic field. From 



_.6S^ = ^ V X [dvA a + l^c^c 



2tt 



(11) 



we obtain J ai = ±e ij djA% = ±e ij djB a = x,y) 
which implies a quantized spin Hall response to the gra- 
dient of external magnetic field Aq = B a with quantized 
spin Hall conductance o s xy = as expected. 

To better understand the topological properties of the 
non-Abelian spin liquid, an effective topological field the- 
ory describing its long-distance and low-energy physics 
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l 



291 would be desired. To do so, we write J a/i 



e ^v\Q^ a a go |_ na £ sa tisfies the continuity condition 
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